AD-A038  595  BALLISTIC  RESEARCH  LABS  ABERDEEN  PROVING  GROUND  MD  F/G  20/4 

AN  INTRINSIC  DESCRIPTION  OF  UNSTEADY  SHOCK  WAVES. (U) 

MAR  77  T W WRIGHT 

UNCLASSIFIED  BRL-1971  NL 


DATE 

FILMED 

' 5 -77 


MICROCOPY  Rl  SOI  111  ION  1 1 SI  ( HARI 

NAIU'NAI  MUM  Ah  * * '• 


REPORT  NO.  1971 


AN  INTRINSIC  DESCRIPTION  OF 
UNSTEADY  SHOCK  WAVES 

Thomas  W.  Wright 


March  1977 


Approved  far  pull  1c  rcliMCi  dlitrlfcwtloc  unlisted. 


If 

Csr 


AO"  DC 

i U 1 APR  26  1977 


UEBEirirr 

A 


USA  ARMAMENT  RESEARCH  AND  DEVELOPMENT  COMMAND 
USA  BALLISTIC  RESEARCH  LABORATORY 

AIN  MEN  ffOVINO  (MOUND,  MARYLAND 


Destroy  this  report  vhen  it  is  po  longer  needed. 
Do  not  return  it  to  the  originator. 


Secondary  distribution  of  this  report  by  originating 
or  sponsoring  activity  is  prohibited. 

Additional  copies  of  this  report  may  be  obtained 
from  the  National  Technical  Information  Service, 

U.S.  Department  of  Conmerce,  Springfield,  Virginia 
22151. 


The  findings  in  this-  report  are  not  to  be  construed 
an  official  Department, of  the  Army  position,  unless 
so  designated  by  other  authorized  documents. 


i 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (mien  Data  Kntered) 

| REPORT  DOCUMENTATION  PAGE  ! m^oUD.SSpletoS"/orii 

p REPQHI-»ty>H«l>  |2.  GOVT  ACCESSION  ~NO.j  3.'  RECIPIENT'S  CATALOG  NUMBER 

BRL  R6p**t~No.-1971 


«.  TITLE  land  lubUtta) 

An  Intrinsic  Description  of  Unsteady  Shock  Waves 

- - S' 


5.  TYPE  OF  REPORT  A PERIOD  COVERED 


6.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHOR!.) 

/Thomas  W./ Wright 
' y /N  / 


6.  CONTRACT  OR  GRANT  NUMBCRfa) 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

US  Army  Ballistic  Research  Laboratory 
Aberdeen  Proving  Ground,  MD  21005 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

US  Army  Materiel  Development  & Readiness  Command 
5001  Eisenhower  Avenue 

Alexandria,  VA  22333  


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  & WORK  UNIT  NUMBERS 


I 


. MONITORING  AGENCY  NAME  » AODRESSf/f  dill, ran!  from  Controlling  Otllco) 

15.  SECURITY  CLASS,  (ol  Wf»  twpvtt) 

Unclassified 

1 5a.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (ol  thl,  Report) 

Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  ebetract  entered  In  Block  20,  ll  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 


This  report  was  prepared  while  the  author  was  the  recipient  of  a Secretary 
of  the  Army's  Research  and  Study  Fellowship  at  Cambridge  University,  Cambridge, 
England,  during  the  academic  year  1974-75. 


19.  KEY  WORDS  (Continue  on  reveree  aide  if  neceeaary  end  Identify  by  block  number) 

Shock  waves,  fluid  dynamics,  growth  equation 


20.  ABSTRACT  (Continue  on  reveree  aide  It  neceeaary  and  Identify  by  block  number) 

In  this  pspcr  the  intent  is  to  obtain,  insofar  as  possible,  an  intrinsic  descrip- 
tion  of  shock  propagation,  that  is  to  say,  a description  such  that  the  motion 
and  the  evolution  of  the  amplitude  of  the  shock  are  determined  by  quantities 
known  on  the  shock  itself.  Such  a description  would  parallel  the  results  of 
geometrical  acoustics,  but  with  significant  differences..  To  begin  with,  a shock 
wave  does  not  generally  lie  on  a characteristic  surface  of  the  governing  equa- 
tions, and  so  it  is  not  possible  to  define  a ray  as  a trajectory  derived  from 


EDITION  OF  1 NOV  66  IS  OBSOLETE 


UNCLASSIFIED  (OVER) 

security  classification  of  this  page  />•<•  EwifWi 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAOEflun  Dili  Bntmnd) 

the  characteristic  equation.  Neither  is  it  possible  to  find  first  the  motion 
and  then  the  amplitude  of  the  shock  since  the  two  are  strongly  coupled. 
Finally,  since  acoustic  disturbances  may  overtake  and  modify  the  shock,  a 
completely  intrinsic  description  seems  to  be  impossible.  Nevertheless,  it 
turns  out  that  a convenient  choice  of  ray  arises  in  a natural  way  from  the 
analysis,  and  that  an  equation  for  the  variation  of  amplitude  along  a ray  can 
be  derived.  Furthermore,  the  coupling  to  the  flow  behind  the  shock  occurs 
only  through  a single  scalar  term,  even  for  the  fully  three-dimensional  case, 
and  the  amplitude  equation  may  be  cast  into  a form  such  that  the  coupling  is 
relatively  weak. 


UNCLASSIFIED 

L 


SECURITY  CLASSIFICATION  of  this  PAGEr»Ti«n  Pur#  Fnfrmd) 


. | ..»IU|||I  LUM.IPU  III. 


TABLE  OF  CONTENTS 


TABLE  OF  CONTENTS  . . 
LIST  OF  ILLUSTRATIONS 

I.  INTRODUCTION 

II.  KINEMATICS 

III.  DYNAMICS 

IV.  TRANSPORT  EQUATIONS  . 

REFERENCES 

APPENDIX 

LIST  OF  SYMBOLS  . . . 
DISTRIBUTION  LIST  . . 


3 

5 

7 

9 

II 

13 

32 

33 
35 
37 


IMWJL 

(TO  *Sfl« 

left  C, 

nmimn  1 




IT 


llHMITtH/MHU»ILin  Won 
, — jai  tl>iL~ul  w 


3 


I . INTRODUCTION 


It  is  often  useful  to  describe  weak  sound  waves  by  the  theory  of 
geometrical  acoustics.  By  this  means,  surfaces  of  constant  phase  in  an 
oscillatory  disturbance  may  be  traced  forward  in  time,  and  amplitudes  of 
the  disturbance,  associated  with  those  surfaces,  may  be  calculated.  The 
heart  of  the  method  lies  in  the  construction  of  a system  of  rays  along 
which  points  of  an  initial  surface  are  mapped  parametrically  by  later  times 
onto  the  successive  members  of  a family  of  surfaces.  Once  the  rays  have 
been  constructed,  amplitudes  are  found  from  a transport  equation,  which 
is  an  ordinary  differential  equation  that  determines  the  variation  along  a 
ray.  The  method  as  described  above  provides  an  asymptotic  solution  in 
linear  acoustics  for  high  frequencies.  By  various  modifications  and 
extensions  the  basic  ideas  may  be  applied  to  pulses,  progressing  waves, 
and  wave  front  expansions.  All  of  this  has  been  conveniently  summarized  by 
Friedlander *,  who  also  gives  several  detailed  examples.  In  nonlinear 
acoustics  acceleration  waves  or  weaker  discontinuities  may  also  be  treated 
geometrically2,  but  here  the  results  are  exact.  In  every  case  mentioned, 
the  ray  equations  and  the  transport  equations  are  only  weakly  coupled  in 
that  amplitude  transport  depends  on  the  rays,  but  the  rays  are  independent 
of  amplitude. 

In  this  paper  the  intent  is  to  obtain,  insofar  as  possible,  an  intrin- 
sic description  of  shock  propagation,  that  is  to  say,  a description  such 
that  the  motion  and  the  evolution  of  the  amplitude  of  the  shock  are  deter- 
mined by  quantities  known  on  the  shock  itself.  Such  a description  would 
parallel  the  results  of  geometrical  acoustics,  but  with  significant  differ- 
ences. To  begin  with,  a shock  wave  does  not  generally  lie  on  a character- 
istic surface  of  the  governing  equations,  and  so  it  is  not  possible  to 

define  a ray  as  a trajectory  derived  from  the  characteristic  equation. 
Neither  is  it  possible  to  find  first  the  motion  and  then  the  amplitude  of 
the  shock  since  the  two  are  strongly  coupled.  Finally,  since  acoustic 
disturbances  may  overtake  and  modify  the  shock,  a completely  intrinsic 
description  seems  to  be  impossible.  Nevertheless,  it  turns  out  that  a con- 
venient choice  of  ray  arises  in  a natural  way  from  the  analysis,  and  that 
an  equation  for  the  variation  of  amplitude  along  a ray  can  be  derived. 

Furthermore,  the  coupling  to  the  flow  behind  the  shock  occurs  only  through 

a single  scalar  term,  even  for  the  fully  three-dimensional  case,  and  the 
amplitude  equation  may  be  cast  into  a form  such  that  the  coupling  is 
relatively  weak. 

In  a previous  paper3  a growth  equation  for  the  shock  amplitude  in  an 
inviscid  simple  elastic  fluid  was  derived.  In  that  paper  only  mechanical 


if.  G.  Friedlander,  Sound  Pulses , Cambridge  University  Press,  Cambridge, 
1958. 

2E.  Varley  and  E.  Cumberbatah,  J.  Inst.  Maths.  Applies,  lj  (1965)  101. 
3P.  J.  Chen  and  T.  W.  Wright,  Meccanioa  (to  appear). 
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effects  were  considered  since  pressure  was  assumed  to  be  a function  of 
specific  volume  alone.  In  addition,  the  fluid  ahead  of  the  shock  was 
assumed  to  be  uniform  and  at  rest.  Finally,  the  analysis  was  based  on  a 
fixed  reference  configuration,  which  is  somewhat  unnatural  for  problems  in 
fluids,  although  the  results  are  perfectly  correct.  This  paper  is  an  out- 
growth of  reference  (3)  and  contains  the  results  of  reference  (3)  as  a 
special  case  with  some  minor  differences  of  notation,  but  there  are  signif- 
icant differences  in  approach.  The  fluid  is  still  assumed  to  be  inviscid, 
but  now  the  pressure  is  taken  to  depend  on  specific  entropy  as  well  as 
specific  volume  although  heat  conduction  is  ignored.  Thus,  the  only  dis- 
sipation effect  originates  in  the  shock  wave  itself.  The  fluid  ahead  of 
the  shock  may  be  in  an  arbitrary  state  with  arbitrary  smooth  motion. 
Finally,  the  analysis  here  does  not  rely  on  a fixed  reference  configura- 
tion, but  rather  makes  use  only  of  the  present  spatial  configuration  at 
the  shock. 

In  the  next  two  sections  the  necessary  preliminaries  are  given.  In 
Section  II,  appropriate  kinematics  to  describe  the  motion  of  a shock  wave 
and  the  evolution  of  shock  amplitude  are  outlined.  The  approach  taken  is 
that  of  the  theory  of  singular  surfaces4,  but  the  idea  of  the  displacement 
derivative  is  not  used.  Rather,  variations  are  sought  along  a convenient 
ray  direction,  chosen  so  as  to  simplify  the  results.  In  Section  III  the 
constitutive  relation,  the  equations  of  motion,  and  the  jump  conditions 
are  reviewed. 

In  Section  IV  the  main  results  are  given.  With  the  derivatives  of 
amplitude  and  entropy  eliminated  in  favor  of  the  derivative  of  shock  speed, 
three  forms  of  the  transport  equation  for  the  variation  of  shock  speed 
along  a ray  are  found.  Each  of  these  forms  exhibits  a different  coupling 
term  with  the  flow  behind  the  shock.  In  each  case  the  exact  value  of  the 
coupling  term  depends  on  the  past  history  of  the  shock  trajectory.  Thus 
the  equations  are  really  functional-differential  equations  in  disguise, 
but  the  explicit  nature  of  the  functionals  is  not  given.  If  it  is  then 
assumed  that  the  values  of  one  coupling  term  are  known  exactly,  then  it  is 
found  that  the  transport  equation  and  a kinematical  relation  between  shock 
speed  and  shock  normal  form  a hyperbolic  system  of  equations.  This  implies 
that  disturbances  to  the  shock  surface  will  propagate  laterally  within  the 
shock  surface  itself.  Finally,  these  results  are  compared  with  Whitham' s 
theory  of  shock  dynamics,  which  was  derived  for  polytropic  gases  by  approx- 
imate means5>6»7>8.  One  form  of  the  equations  given  here,  considered  as  a 

A.  Truesdell  and  R.  A.  Toupin,  "The  Classical  Field  Theories,"  Flugge's 
Handbuoh  der  Physik,  III /l.  Springer,  Berlin-Gottingen-Heidelberg , 1960. 

$G.  B.  Whitham,  J.  Fluid  Mech.  2j  (1957),  145. 

6G.  B.  Whitham,  J.  Fluid  Mech.  5,  (1959),  369. 

7G.  B.  Whitham,  J.  Fluid  Mech.  3lj  (1966),  449. 

6G.  B.  Whitham,  "Linear  and  Nonlinear  Waves,"  John  Wiley  and  Sons  (Wiley- 
Interscience) , New  York  - London  - Sydney  - Toronto,  1974. 
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hyperbolic  system  and  specialized  to  polytropic  gases,  gives  the  same 
speed -of^propagatdron-wirthin  the  surface  as -Whit-ham’ s theory  -for  ail  Mach 
numbers. 


II.  KINEMATICS 

A shock  wave  may  be  described  as  a propagating  singular  surface  across 
which  the  pressure,  density,  entropy,  and  particle  velocity  are  all  dis- 
continuous. There  are  several  convenient  mathematical  representations  of 
the  surface.  For  example. 


£(x,  t)  = o 

t - t (x)  =o  (2 . 1) 

1 2 . 

X = x(u  , u , t) 

In  these  equations  x is  the  ordinary  Euclidean  position  vector  in  a fixed 

1 2 

Cartesian  coordinate  system,  t is  time  and  u , u are  Gaussian  surface 
coordinates.  The  normal  in  the  direction  of  propagation,  n,  and  the 
normal  velocity,  v,  may  be  computed  from  any  of  the  three  representations. 
The  first  two  give 


Z : I = t : 
>X  ,t  ,X 


-1  = o:  -v 


(2.2) 


where  the  comma  denotes  partial  differentiation.  With  the  surface  coordi- 
nates orientated  correctly  for  sign,  the  third  representation  gives 


n = x xj2  / x xj2| 
v = X . n 

9 L 


(2.3) 


Equation  (2.3)~  may  be  derived  by  inserting  (2.1)  into  (2.1)  , differ- 

Z 1 2 ^ 
entiating  with  respect  to  t holding  u , u fixed,  and  applying  (2.2). 


The  function  x may  be  regarded  as  a coordinate  transformation  from 
12  12  3 

u,u,ttox,x,x.  Its  partial  derivatives  have  special  significance. 
Basis  vectors  on  the  surface  are  given  by 


»r  * xir  ; r ■ l,  2 


(2.4) 


Greek  capital  letters  will  always  be  used  for  components  of  surface  vectors 
and  tensors,  and  therefore  will  always  have  the  range  1,  2 unless  it  is 
stated  otherwise.  The  derivative  with  respect  to  time  will  be  called  the 
propagation  vector. 


b = 


(2.5) 
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Clearly,  b need  not  be  normal  to  Ap  and  A2  since  infinitely  many  representa 

— tions  of -the  ^orm  (2rl)~could  be  made.  Tn~fact  , orre  t>f  the  principal  tasks 

of  this  paper  is  to  choose  b in  a convenient  and  simplifying  manner. 

The  surface  metric  tensor  and  reciprocal  basis  follow  in  the  usual 
way  from  (2.4) 


afa  " Ar 


aa 


A = det(ApA)  (2.6) 

APA  = (A^)'1 
, r . ta  . 

A = A A. 

'A 

The  coordinate  transformation  will  be  invertible  if  its  Jacobian  determi- 
nant, J,  does  not  vanish. 


(Ap  x A2)  • b = |det  ApA| 2 n ■ b 


A 2 v ^ o 


(2.7) 


It  is  now  easily  checked  that  the  inverse  of  the  Jacobian  matrix  is  the 
following 

aur  /,  1 J 

3 T = (1  - - D » b^)  A 


(2.8) 


3u" 

Sx 


0 

v 


In  (2. 8)_  the  substitution  u = t has  been  made  to  indicate  that  time  is 

^ 12 
to  be  paired  with  the  surface  variables  u , u and  not  with  the  spatial 

12  3 

variables  x , x , x . For  clarity  this  practice  will  be  followed  in  the 
sequel,  and  the  symbol  t will  only  stand  for  time  in  conjunction  with 
12  3 


If  a function  is  defined  only  on  the  shock  surface,  it  is  most  natur- 

12  3 

ally  regarded  as  a function  of  u , u , u . Yet  with  the  aid  of  (2.8)  it 
still  makes  sense  to  compute  the  spatial  gradient  of  such  a function. 

12  3 

Thus,  if  the  amplitude  of  the  shock  wave  is  a = a (u  , u , u ) , its  gradient 
is  given  by 


9a  /.  1 u\*T  n 

7-  = a.  1 - — n a b)  A + aT"*- 
3x  r - V ' ~ ~ ,3  v 


(2.9) 


L 
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and  the  directional  derivative  along  b is  given  by 


b.|=a)3  (2.10) 

The  notation  a^,  indicates  covarient  surface  differentiation.  Other 

functions,  such  as  the  pressure,  are  defined  for  all  x,  t,  but  the  limit- 
ing values  on  the  shock  surface  are  of  particular  interest.  For  such  a 
function,  f(x,  t) , let  f and  f"  indicate  the  limits  taken  from  the  region 
ahead  of  or  behind  the  shock  respectively.  Thus,  we  have 


f*  = lim  f(y,  t(x)  * e)  (2.11) 

e+o 

The  difference  or  jump  between  the  two  sides  is  indicated  as  follows. 


[f]  = f"  -f+ 


(2.12) 


In  taking  the  gradient  of  such  a function;  one  must  be  clear  as  to  the 
order  of  the  differentiation  and  the  limit  on  the  surface.  From  (2.11) 
and  (2.2),  it  follows  that 


(f )*  . = (f.)1  ♦ — (f  J- 
,1  K 1 v , t ' 

n . 

if]  . = [f  . ] + — [f  j 

1 J ,1  ,1  V ,t 


(2.13) 


Formulas  (2.9)  and  (2.10)  apply  to  the  functions  f and  f as  well. 


III.  DYNAMICS 

In  this  section  one  version  of  the  standard  equations  of  inviscid 
fluid  mechanics  is  reviewed.  Let  e,  the  specific  internal  energy  of  the 
fluid,  be  regarded  as  a function  of  specific  volume,  u,  and  specific 
entropy,  n.  We  have 


e = e(u , n)  (3.1) 

In  terms  of  this  thermodynamic  potential  the  pressure,  p,  and  temperature, 
6,  are  given  as  follows. 


0 = e 

n 

The  subscripts  denote  partial  differentiation. 

In  regions  adjacent  to  a shock  wave,  the  conservation  laws  for  mass, 
linear  momentum,  and  energy  may  be  expressed  by  partial  differential 
equat ions. 


u = u u. 

1,1 


-p  . = pf.  + pu. 
r,l  1 ^ 1 


(3.3) 


n « 0 

The  third  equation  is  a reduced  form  for  the  conservation  of  energy  in 

the  absence  of  heat  conduction  and  supply  and  makes  unnecessary  a separate 

inequality  for  entropy.  The  particle  velocity  is  u,  the  density  is  p,  and 

the  body  force  is  f.  Latin  indices  are  associated  with  Cartesian  vectors 

and  tensors.  The  dot  denotes  the  material  time  derivative,  e.g. 

u = u ^ + u . u. . 

,t  ,1  1 

Across  a shock  wave,  the  conservation  laws  must  be  expressed  by  jump 
conditions.  It  is  useful  first  to  define  , the  normal  shock  speed 
relative  to  the  fluid  on  either  side  of  the  shock. 


IT  = v - u.  n.  (3.4) 

The  conservation  laws  for  mass,  linear  momentum,  energy  in  a reduced  form 
and  entropy  may  now  be  stated  as  follows. 


[pU]  = o 

[p]n  = p+  U+  [u. ] 

+ + + (3-5) 

P U [e]  = !*(P  + P ) [u.]  nA 

[n  ] > o 

It  is  assumed  that  p^  is  always  negative,  but  p^y  need  not  be  positive  for 

the  development  of  the  theory.  It  is  assumed  that  all  shock  waves  to  be 
examined  not  only  obey  equations  (3.5),  but  are  stable  as  well.  The 
questions  of  necessary  or  sufficient  conditions  for  stability  will  not  be 
considered . 

In  a shock  wave  U f o.  Thus  equation  (3.5)  shows  that  the  jump  in 
particle  velocity  must  be  normal  to  the  wave,  and  therefore,  it  may  be 
expressed  as  follows. 


[u. ] = a U+  n. 
i 1 l 

- + + 
oru.  =u.  + a U n. 

li  l 


(3.6) 


U has  been  introduced  as  a convenient  normalizing  factor,  and  a is  called 
the  amplitude  of  the  wave.  It  now  follows  directly  from  (3.4),  (3.5)  , and 
(3.6)  that  1 


u = u+  (1-a) 
U"  = U+  (1-a) 


(3.7) 
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(3.8) 


and  from  (3.5),,  and  (3.6)  that 

p*  (uV  - Jfi 

From  (3.7)^  it  is  clear  that  the  amplitude,  a,  is  a measure  of  the 
compression  (or  dilation)  produced  by  the  shock. 


The  entropy  behind  the  shock  can  be  written  in  terms  of  b,  defined  as 
follows 

n"  = n+(l+b)  (3.9) 

Of  course  it  is  required  that  b > o.  Since  e and  p were  originally  taken 
to  be  functions  of  u and  o,  it  is  clear  that  immediately  behind  the  shock 
they  may  now  be  regarded  as  function  of  a and  b with  u+  and  n + as  (known) 
parametric  functions.  Equation  (3.8)  is  one  relation  among  the  variables 
(a,  b,  U+)  and  (3.5),  with  (3.6)  provides  another. 

P+  [e]  = *2(p+  + p")a  (3.10) 


In  some  cases  it  may  be  convenient  to  use  (3.8)  and  (3.10)  so  as  to  elimi- 
nate two  of  the  variables  (a,  b,  U+)  in  terms  of  the  third.  For  example. 


The 

+ 

u . 
b = 


da 


<ap  -2  p+  e > 

— It-: -Hr  u+  du+ 

en'  Po' + ^u+n 


(3.11) 


db  = dU+ 

n e 

n 

solution  to  these  differential  equations  will  depend  parametrically  on 
nf,  and  leads  to  the  well  known  results  that  for  small  amplitudes 

0(a3). 


IV.  TRANSPORT  EQUATIONS 

If  equations  (3.8)  and  (3.10)  could  be  supplemented  by  a third 
independent  relation  among  the  three  variables  a,  b,  U+  with  time  as  a 
parameter,  then  the  evolution  of  these  quantities  during  shock  propagation 
could  be  completely  determined,  at  least  in  principle.  If  such  a relation 
is  available,  it  must  exhibit  dynamical  coupling  with  the  flow  behind  the 
shock  since  none  exists  in  either  (3.8)  or  (3.10),  and  yet  such  coupling 
must  be  present. 

A shock  wave  connects  two  adjacent  regions  of  smooth  flow  where 
equations  (3.3)  hold.  An  attempt  to  obtain  the  desired  third  relation  can 
be  made  by  manipulating  the  jump  in  (3.3), 


A 
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-pu'  (u,ir  -Pn"  ^fiJ'  +pu  (u,i)+  +pn+  ^tiy 
+(p  -p  J f.  = p u.  - p u. 


(4.1) 


The  procedure  now  followed  is  to  express  as  many  quantities  in  (4.1)  in 

terms  of  a,  b,  U and  u , n , u.  as  possible.  The  terms  p . p and  p4 , 

Pn  already  have  the  required  dependence  through  (3.7)  and  (3.9).  The 

process  of  converting  (u  . ) and  (n  .)  into  the  required  form  is  lengthy 

but  straightforward,  and  requires  little  more  than  repeated  application 
of  (2.13),  (3.3)j,  and  (3.3)^  to  (3.7)  and  (3.9).  Begin  with  (2.13) 

applied  to  specific  volume. 


+ + “i  + 

(u  . )~  = (u“)  . - — (u  )- 
,1  ,1  V v ,t' 


(4.2) 


Either  the  top  sign  or  the  bottom  sign  is  to  be  chosen  consistently.  From 
the  definition  of  material  derivative  we  also  have 


. + • + 4-  + 

(u  ) = (u)  - (u  .)  u. 

> L ,11 


Equations  (4.2)  and  (4.3)  combine  to  give 


{ 1 +\  + -4-  ni  . + 

{6ij  - v ni  uj"J  (0,jr  = 4 - — (tf 


,1  V 

The  matrix  on  the  left  hand  side  is  easily  inverted. 


(u  J = (V;  + 


1 


j.  . ■ — - n,  u . 
ki  ± k i 


n . ) 

— I 


(u1)  . - — (0) 

,i  v \ 


(4.3) 


(4.4) 


(4-5) 


For  the  minus  signs  in  (4.5)  and  with  the  use  of  (3.3).,  (3.6)2>  (3.7)  , 
and  (3. 7) ^ this  becomes 

(u  k)_  = i6ki  + ~ (ui  + a U+ni)| 

’k  ( kl  U+  (l-a)  1 1 j 


n . 

(u+)  . (l-a)-0+a  . - -i  u+ (1 -a)  (u 

>-*-  J 9 J 


(4.b) 


The  last  term  in  (4.6)  may  be  expanded  with  the  aid  of  (2.13).  and  (3.6)~ 
to  give  (4.7) 


n . 

l 


(uj.jr  “ (“j  * 9 U*  V.j  - 4T  (uj,tr  • 


(4.7) 


and  the  last  term  in  (4.7)  may  be  expanded  to  give  (4.8). 

U*(I-a)  _ 

v ')■ 


(u.  *)  ~ Vr  = (*i ) -(Uj1  + a U n^ ) ■ 


{‘"j*1.!  * 19  "*  nj>,l} 


(4.8) 
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To  arrive  at  (4.8),  the  definition  of  material  time  derivative  as  well  as 
(2.13).,  (3.4),  and  (3.6)2  have  been  used.  Inspection  of  (4.6)  - (4.8) 

shows  that  spatial  derivatives  must  be  taken  on  two  types  of  quantities: 
those  defined  only  on  the  surface,  and  those  defined  in  the  region  ahead 
of  the  shock  but  restricted  to  the  shock  as  in  (2.11).  Accordingly,  the 
derivatives  must  be  interpreted  by  (2.9)  or  (2.13). 

The  term  (u+)  . may  be  found  from  (4.4)  with  (3.3)  substituted. 

+ 1 

n . u . n . 

(u+)  . = (6.  — — 2-0  (u  .)  + ♦ — u+  (u.  .)+  (4.9) 

,1  ij  v ,r  v j,y 

Similarly,  for  (u.+)  . we  have 

1 ’ n.  u,  ) n. 


* = \6^  - ! K J + + 


(4.10) 


Finally,  with  (4.7)  - (4.10)  equation  (4.6)  reduces  to  the  following 


<u,i>'  ■ ! 


u n.n . 


(l-a)6  + a n.  n l (u  )+  + - -j-* [u  ] 

1 J J f ,J  (U ) (1-a)  J 


!.  n . u . 

<5 . . + r~-  n . n . + 2 — — 2 

X-a  1 J U+ (1-a) 


u+  n.  ( U+  n,  + u+,  1 

\ k k J (U  + ) , - n.  u+  a n.  . 

U+  t U+  (1-a)  ) ’k  1 


(4.11) 


- — — — n.n.n,  (u.  , ) + 

U+(l-a)  1 J k 

r 

The  term  -n . . in  (4.11)  is  equal  to  Br , which  is  twice  the  mean  surface 
J.J  r 

curvature.  That  this  is  true  may  be  seen  from  (2.9)  applied  to  n and  the 

r r 

formulas  n • A = -B„  and  n‘Dr  = n • D , = o.  The  fact  that  n = o 

,r  r , 

simplifies  the  reduction  of  (n  .)  , but  the  process  is  exactly  the  same. 
The  result  is 


(n  •)'  = (1+b)  Is-.  + r^-n.n.l  (n  -)  + 
n,i  \ ij  1-a  1 jf  U,]J 

( n . u . 1 

+n  <6.  . + n.n.  + — — -2 — \b  . 

\ 1J  1‘a  1 J l)+(l  -a) ( 

The  body  force  and  acceleration  terms  may  be  rewritten  as 

[pU  - 1 - [p]  f-  = [u . ] + (U.+  - f.) 

1 1 lt  1 i 1 -a  1 1 1 1-a  1 1 

■ hf  iai>  - 


(4.12) 


(4.13) 
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Equations  (4.11)  - (4.13)  should  be  inserted  in  (4.1)  to  complete  the 
desired  reduction. 

Rather  than  writing  out  the  whole  equation,  it  is  more  convenient  to 
examine  the  normal  and  surface  components  separately.  But  even  before 
doing  that,  note  that  a certain  combination  of  derivatives  of  the  amplitude 
arises  when  taking  the  normal  component  of  (4.11).  That  is 

(U+  n.  + u.+)  a . 

1 i l ’ ,i 

The  same  combination  occurs  for  U+  and  b.  No  choice  has  yet  been  made  for 
the  propagation  vector  in  (2.5),  so  we  now  set 


t?  = x =U+n+u+  = U n+u  (4.14) 

9 ^ 

where  the  last  equality  follows  from  (3.6)2  and  (3.7).,.  From  (3.4)  it  is 

clear  that  this  choice  satisfies  (2.3)2-  In  effect,  this  one  feature  of  the 

surface  parametrization  (2.1)^  has  been  chosen  strictly  for  convenience. 

Any  other  choice  that  satisfies  (2.3).,  could  also  be  made  to  serve.  Equa- 

*- 

tion  (4.14)  defines  a family  of  rays  along  which  the  surface  may  be  thought 
to  propagate.  Now  with  the  aid  of  (2.10)  and  (4.14),  the  normal  component 
of  (4.1)  may  be  set  down.  The  plus  sign  is  always  understood  to  apply  to 
all  terms  unless  the  minus  sign  is  explicitly  indicated.  This  convention 
will  be  adopted  in  the  sequel. 


2 op  a _ 

u ,3 


+ 


aup 

u 

U 


U 


,3 


. - ..  /uPu  u\  , 

= -a(l-a)oPu  Un.f.  n.[u.] 

(4.15) 

+U  {Po"  (1'a)  -pu}  <%>  + lJ{Pn’  (1+b)‘Pn}(n,n) 

-aup  (n.u.  ) 

u j J ,n 

The  subscript  n denotes  the  normal  derivative,  e.g.,  (u  ) = n.(o  .). 

Two  other  equations  for  the  3-derivatives  of  a,  b,  U may  be  found  by 
differentiating  (3.8)  and  (3.10)  and  applying  the  following  two  simple 
ident it ies . 


(o)  = U(o  n)  * u (u  ) 

9 J 9 11  J 9 J 

<">.3  = U<\n> 


(4.16) 


The  first  of  these  is  derived  from  (2.10),  (4.14)  and  (4.9),  and  the  second 
is  derived  similarly.  The  two  new  equations  may  now  be  written  down. 


1(> 


-(°PU  + pU  )a,3  + °Pn  b,3  ' 2apUU,3 

= - {(l-a)pu“  - Py  ♦ a(PU)2}  {llu^  ♦ uUjJ| 

- {(l+b)pn-  -Pn}U^,n 

- J {upo~  + plj2}  a, 3 +{l3nPn'  " PnCn’}  b,3 

= '{I  [Pu  + Pu’  (1-a)]  + a(14)  (pU)2}{Uu,n  +uuj,j} 
' {I  [pn  + Pn"  (1+b)]  + p[en  " en"  (1+bJ}  Ur>,n 


(4.17) 


Rather  than  solving  directly  for  a ,,  b , and  U , from  (4.15)  and 

, O , J f -J 

(4.17),  it  is  convenient  to  regard  a and  b as  functions  of  U,  as  in  (3.11), 
Then  (4.17)  may  be  solved  for  a , and  b - in  terms  of  U . 

f O , «>  > *) 


,3 


[l  - If1!!  . 

1-M  I ’Pe"  j 


en  Pu  (1'p  5 


(2  T1  " lz) 


P„  (1-v  ) 


ll 


(4.18) 


,3 


a2U 

ne. 


U 3 + 

t 


TT  {Ih  - "2} 

ne  ’ 


n n 

In  (4.18),  1^  and  I are  the  right  hand  sides  of  (4.17)  and  (4.17)2 

respectively.  Since  1^  and  vanish  for  homogeneous  conditions  ahead  of 

the  shock,  they  will  be  called  the  inhomogeneous  terms.  In  addition,  p, 
the  Mach  number  behind  the  shock,  has  been  introduced.  We  have 

2 2 
= . (pu) 


2 

U = 


7) 


and 


(4.19) 


■f  j • ■ •f 


for  the  Mach  number  ahead  of  the  shock.  Here  c and  c are  the  acoustic 
speeds  relative  to  the  flow  ahead  of  and  behind  the  shock,  respectively. 
1 + 2 — 

We  have  (p*c_)  = -p()  . 

17 


With  the  aid  of  (4.18),  (4.15)  becomes 


(4 • 20) 

= -a(l-a)  ^n.  . + 0 • [u]  + I 

U ’ V 

I is  the  net  inhomogeneous  term  from  (4.15)  and  (4.18)  and  is  given  in 
full  in  the  Appendix.  Equation  (4.20)  together  with  (4.14),  repeated  here, 

x,  = Un  + u (4 . 14) 

describe  the  instantaneous  motion  of  the  surface.  At  a given  instant  of 
time,  if  the  location  of  the  surface  is  known,  then  the  normal  q and  mean 

r 

curvature  B^,  are  known.  Furthermore,  if  the  distribution  of  1)  over  the 

surface  is  known,  then  the  distributions  of  a and  b are  known.  Finally, 
since  u,  n,  and  u are  known  ahead  of  the  shock,  the  only  remaining  term  in 
(4.14)  and  (4.20)  that  is  not  known  is  the  single  scalar  quantity  n • [y]. 
This  final  quantity  provides  the  only  coupling  of  the  shock  with  the  flow 
behind  the  shock  in  spite  of  the  fully  three-dimensional  nature  of  that  flow. 

Equation  (4.20)  shows  that,  aside  from  inhomogeneities,  it  is  the 
competition  between  the  mean  curvature  and  the  coupling  to  the  rearward  flow 
that  determines  whether  U is  increasing  or  decreasing  with  time,  and  hence, 
through  (4.18)  whether  the  amplitude  is  increasing  or  decreasing. 

In  a similar  way,  with  (4.11)  - (4.13)  inserted  in  (4.1),  equations 
connecting  the  transverse  derivatives  of  a,  b,  and  l)  may  be  found  by 
resolving  (4.1)  along  the  base  vectors  in  the  shock  surface. 

upo'  a|r  ' npn’  bJ  r = T^a  *r  * lpJ 

+ {cl'a)  pu“  ' rbp4  hr3  t4,21) 

+ {n+b)pn‘  - T^pn}  (V 

With  a and  b determined  over  the  whole  surface  from  the  solution  of  (4.14) 
and  (4.20),  this  equation  determines  the  transverse  components  of  acceler- 
ation behind  the  shock.  It  may  be  thought  of  as  a compatibility  condition 
for  it  indicates  that  acceleration  may  not  be  arbitrarily  prescribed  behind 
a shock  wave. 


l + 3y 


2 M 2. 

u (l-y  ) 


a - 


IV 

iV 


pe 
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Other  Forms  of  the  Transport  Equation 


Although  (4.20)  is  exact  as  it  stands,  it  may  be  cast  into  other  forms 
by  choosing  other  terms  to  provide  coupling  with  the  rearward  flow.  A 
discussion  of  the  growth  and  decay  of  plane  shocks  in  an  elastic  solid9 
used  the  strain  gradient  as  the  coupling  parameter;  and  a similar  discussion 
concerning  curved  shocks  in  an  elastic  fluid  (see  reference  3)  used  the 
totally  normal  component  of  the  second  gradient  of  deformation,  i.e.,  the 
term  c = n.N  N_(x.  where  the  deformation  is  given  as  a function  of 

material  coordinates  and  time,  x = xx(Xa,t),  and  the  spatial  and  material 
normals  are  n and  N respectively.  In  the  present  context,  the  quantity 
most  closely  related  to  these  two  is  the  normal  density  gradient  behind 
the  shock  or  alternatively  the  normal  derivative  of  specific  volume  behind 
the  shock.  The  following  derivation  uses  only  kinematic  identities  and 
jump  conditions  to  show  the  relationship  between  the  normal  component  of 
acceleration  and  the  normal  derivative  of  specific  volume  behind  the  shock. 


Begin  with  the  following  identity,  which  uses  (3.6)2  and  (4.14). 


nj  (U‘  nk  + uk"3  (uj"\k 

= n5  (U+  nk  + uk+)  (u.+  ♦ a U+  n.)>k 


(4.22) 


With  the  aid  of  (2.13)^  and  (2.10),  the  terms  of  (4.22)  may  be  expanded  and 
regrouped . 


n.[i.]  ♦ njnkU_ [Uj  k] 

= a U+  n.nk(uj>k)+  ♦ (a  U+)>3 

Next,  note  that 


(4.23) 


n.n.(u.  .)  = (u.  .)  - (6--  -n.n.)  (u.  .) 

1 1 i,J  1,1  11  1 J i,J 


(4.24) 


The  last  term  in  (4.24)  involves  only  surface  derivatives.  Since  the 
transverse  component  of  velocity  is  continuous  across  a shock,  we  have 

Ta^  u.l  = o = aF  [u.  ] + Ar  A.  [u.  .] 

L 1 ij  ,r  i,r  1 1 jr  i,r 

which  may  be  rewritten  as 

-a  U+  n.  . + (6-.  - n.n.)  fu.  .]  = o (4.25) 

1,1  voij  1 y 1 i,jJ 

From  conservation  of  mass  (3.3)^  and  a rearrangement  of  terms  we  have 
(ui,i)"  = p’  (p)’  = p'  {(u)",3  • U’  ni(u,ir} 


J.  Chen  and  M.  E,  Gurtin,  Int.  J.  Solids  Struc. 


(1971),  5. 
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which  may  be  rewritten  as 


- -nr*  <Vi>'  - i»u>'  i».»i 


(4.2b) 


where  (3.3)^,  (3.7)^  and  (4.1b)j  have  been  used,  liquations  (4.23)  - (4.2b) 
may  be  combined  to  give 


nJu  ] = (l-a)plT  [u  J + 2 U a 

J J P 11  » *5 

+ a U + a(l-a)  U2  n . 
» ^ ^ * 1 

+ a n.n.  U u.  . 
j k j ,k 


(4.27) 


where  all  terms  on  the  right  hand  side  carry  the  plus  sign.  With  the 
substitution  of  (4.27)  and  use  of  (4.18)j,  equation  (4.20)  becomes 


1+3  u 
l-u2 


1-3U2  Pn 


l-u 


pe. 


,3 


(4.28) 


= -a(l-a)u2  n.  . + — -y-  (l-a)piniu  ] + I 
i,i  2 1 ,n 


where  I stands  for  all  left  over  inhomogeneous  terms  and  is  given  in  full 
in  the  Appendix. 


Another  possible  coupling  term,  suggested  by  Whitham's  papers  on  shock 

dynamics  (see  References  5,  6,  and  7),  is  the  combination  (p  , ♦ pcn.u.  ,)’. 

» t 1 1 t t 

The  derivative  with  respect  to  t'  signifies  the  time  variation  as  seen  by  an 
observer  who  is  instantaneously  at  rest  with  respect  to  the  flow  ahead  of 

the  shock.  That  is  to  say,  + u.  — — . The  following  identities 

dL  dt  1 3X  ^ 

will  prove  useful. 


(f  t,r  = (f)-  - [u  ] (f  )■ 

(f)  = (ll+  n + u +)  (f") 

t O XI  ,1 

= (f  t,)‘  + U+  n (f  )• 

yL  1,1 


(4.29) 


It  is  possible  to  work  out  a relationship  between  n.(u^)  and 

(p  , + pcn.u.  ,)  for  substitution  in  (4.20),  but  it  is  simpler  to 
, t l l , t 

proceed  directly  as  follows.  With  the  aid  of  (4.29)j  we  have 

(P,t»  + pCniUi,t'r  = Pu"(u  Ui,ir  - tuiJ(P,ir 


(PC)*  ni  (-  P\i  + fi) 


(4.30) 


- (pc)'  n I u . | (u  )' 
1 J 1 > J 


JO 


where  (3.3)  has  been  used  to  eliminate  p and  u..  The  acoustic  speed 

- - 2 1 

behind  the  shock  is  given  by  (p  c ) = -p  so  that  rearrangement  of  terms 

in  (4.30)  gives 

(P.f  + pcniui,tf>“ 

= -(c"n.  + [u.j)(p(i)"  + (pc)"  n.f. 

(4.31) 

- (pc)  ni (c  n^  + [u.])  (Ui  j) 

- (pc2)  (6 . . - n.n . ) (u.  . ) 

ij  i J i,) 

fince  [u]  = aU+n,  the  first  and  third  terms  involve  only  normal  derivatives, 
which  may  be  reexpressed  by  (4.29)2»  and  the  last  term  may  be  rewritten 

with  (4.25).  With  rearrangement  of  terms  and  use  of  (3.7),  equation  (4.31) 
becomes 


(p,t-  + ^Vi.t'1' 


* {(p_),3  + (pc)"  V'O.s} 

+(pc)  n.f.  - (pc2)  |(5. .-n.n.)  (u.  . ) + + aU+  n.  .1 

li  v \ ij  ij  i,j  1,1  I 


(4.32) 


When  this  equation  is  multiplied  through  by  (l-a)U  /c  , and  (p)  , and 

9 ^ 

(u^  ) j are  expanded  with  the  use  of  (3.6)2,  (3.7)j,  and  (3.9),  we  have 


[l-a(l-u)] j(-uPu"  ♦ — ) a 3 + npn"  b>3  ♦ — U>3 


-a(l-a)  £ — n + (l-a)(l-y)(p  , + pcn.u.  , 


(4.33) 


-li-ad-u)]  j(l-a)Po-  u>3  ♦ (l-b)?^  n>3  +£_n.(u.)j3j 
„ , (1-a)  U2 

♦pUn.f.  - - (6-.  - n.n.)u.  . 

M 1 1 u ij  1 J i.J 

When  combined  with  (4.18),  equation  (4.53)  may  be  reduced  to 


j{^  j - »2Ju,3 


(4.34) 


= -a(l-a)  n • 4 (P>t . + Pcni«i,t.>_  + 1 


where,  as  before,  I stands  for  the  inhomogeneous  terms  and  is  given  in 
full  in  the  Appendix. 
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(4.35) 


Equations  (4.20),  (4.28),  and  (4.34)  all  have  the  form 

A U _ + B U2  n.  r = C + I 

,3  l i,r 

where  A and  B are  nondimensional  coefficients,  C is  a term  that  couples 
the  shock  to  the  rearward  flow,  and  I represents  the  remaining  inhomogeneous 
terms.  The  terms  A,  B,  C,  and  I have  different  forms  in  each  version  of 
(4.35),  of  course.  In  every  c.ise  A,  B,  and  I are  known  functions  of  U,  n, 
x,  and  t.  The  coupling  term  C may  be  further  decomposed  into  a product  of 
which  the  first  term  is  also  a known  function,  but  the  second  term  contains 
derivatives  of  flow  quantities  p , u , or  u . These  derivatives  are  not 
known  in  terms  of  quantities  intrinsic  to  the  shock  wave,  but,  rather, 
depend  on  the  history  of  the  flow.  In  particular,  they  depend  on  initial 
and  boundary  data  where  the  shock  wave  itself  is  one  boundary,  although 
there  may  be  others.  Thus,  (4.20),  (4.28)  or  (4.34)  should  more  properly 
be  regarded  as  functional-differential  equations  rather  than  differential 
equations . 


Nevertheless,  let  us  assume  for  the  moment  that  one  of  the  coupling 
terms  is  known  exactly.  Then  (4.35)  and  the  equation 


n . 

i. 


, + Ar  U 

„■>  l 


,r 


-n . 
J 


) 


+ 


r 

A. 

l 


(4.36) 


form  a hyperbolic  system  in  the  unknowns  U,  n.  Equation  (4.36)  is  a purely 
kinematic  relation.  It  has  been  obtained  by  noting  that  since  n • A = o, 

we  have  • Ar  = -n  • A^.  But  A^.  = g>r3  = x>3r  = b>r  and  n ' 0 3 

so  (4.36)  follows  immediately  with  the  use  of  (4.14). 


1 


3. 


The  function  <t>(u",  u“,  u*')  = o for  the  characteristic  surfaces  of 
(4.35)  and  (4.36)  must  satisfy  the  equation 


./{ 


2 TA 
A($  ) - BA 

* ** 


,r*,4  = 0 


(4.37) 


where  A represents  the  contravarient  components  of  the  metric  tensor  for 

the  shock  surface.  The  case  <p  „ - o is  spurious,  for  it  corresponds  to  the 

» ^ 

requirement  that  surface  derivatives  of  U and  n are  continuous  on  the  shock 
surface.  Hence,  the  remaining  term  in  brackets  in  (4.37)  must  vanish, 
which  corresponds  to  wave  propagation  within  the  shock  surface  itself. 

That  is  to  say,  disturbances  to  the  shock  surface  will  spread  as  a wave 
within  the  surface  itself.  Equation  (4.37)  indicates  that  the  speed  of 
propagation  is  independent  of  direction.  If  s is  the  speed  and  1 is  any 
unit  vector  that  is  tangent  to  the  shock  surface,  we  have 


<t> 


,3 


Ai  + ,r  = 5 


l . 


(4.38) 
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Strictly  speaking,  equations  (4.35)  and  (4.36)  should  be  supplemented  by 

„ 12  3 12  3 

other  equations  since  the  functions  x(u  , u , u ) and  Ap(u  , u , u ) are 

not  known  before  hand,  yet  they  are  present  in  the  coefficients  and  the 

right,  hand  sides  of  (4.35)  and  (4.36).  Appropriate  supplementary  equations 

are  x „ = b and  A _ = b . Now  the  list  of  unknown  functions  to  be  deter- 
',3  ' ' T, 3 -,r 


mined  is  U,  n,  jc,  and  A 


If  C is  again  regarded  as  known,  the  full  set 


of  equations  is  hyperbolic  as  before,  and  the  only  nonvacuous  characteristic 
condition  is  still  that  the  bracketed  terms  in  (4.37)  must  vanish.  Of 
course,  it  should  be  shown  that  a solution  to  the  full  set  of  equations  is 
self-consistent  since  U,  n and  A must  all  be  derivable  from  x,  but  that 
proof  will  not  be  attempted  here! 

For  arbitrary  levels  of  shock  amplitude  the  speeds  calculated  from 
(4.20),  (4.28),  or  (4.34)  will  not  be  the  same  since  the  functions  A and  B 
are  not  the  same  in  the  three  cases.  However,  for  weak  shocks  the  three 
calculated  speeds  are  nearly  the  same  and  all  vanish  for  infinitesimal 
shock  amplitudes.  Since  y^l  and  a-*o  for  weak  shocks,  the  speed  tends  to 


(4.39) 


Entropy  changes  may  be  ignored  for  weak  shocks  so  (4.39),  expressed  as  a 
function  of  amplitude  becomes 


r + + 


HPo  I 

If  p ^<o,  then  a<o  and  (4.40)  must  be  modified  accordingly. 


(4.40) 


For  shocks  of  stronger  amplitude  it  is  not  easy  in  general  to  compare 
the  wave  speeds  s for  the  three  cases  since  the  internal  energy  function 
and  its  derivatives  have  not  been  specified.  The  polytropic  gas  is  an 
important  special  case  for  which  an  explicit  comparison  is  easily  made, 
however.  The  specific  internal  energy  in  this  case  is  given  as  follows. 


Y-l 


exp 


n-n. 


(4.41) 


where  e , u , and  n are  reference  quantities,  c is  the  specific  heat  at 
constant  volume,  and  y is  the  ratio  of  specific  heats  (c  , the  specific  heat 
at  constant  pressure,  divided  by  c ).  Formulas  for  pressure,  temperature, 
or  other  derivatives  of  e are  easily  found. 
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p " -eu  ‘ V * 
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(4.42) 
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y-i 


= £_ 
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Y , . 2 

Pu  = ' u P = - (PC) 

By  manipulating  the  jump  conditions  (3.8)  and  (3.10)  in  conjunction  with 
(3 . 7 ) ^ and  (4.42)^,  one  may  easily  find  formulas  for  p~  and  a in  terms  of 

the  Mach  number  m. 


a = 


P = 


easily  calculated. 
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(4.43) 


Pn’/Pen'  are 


(4.44) 


The  relations  among  m,  p,  and  a are  shown  in  Figure  1 for  Y = 1.4.  With 
these  formulas  the  speed  of  disturbances  on  the  shock  surface  always  has 
the  form 


(4.45) 


where  the 

wave  factor 

(4.20): 

A,  (m)  = 2 — 
1 Y 

(4.28): 

A2(m)  = ij  ' 

M 

(4.34); 

A,(m)  = 

-12  y-3  1 

' 3 ‘ ^ 7 

|l  * 3i.2  . 2(1-30  7)  Tri 

[ u 1 


(4.46) 
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2(l-Q(l-u) 
2+p2  (y-1)  . 


These  functions  are  plotted  in  Figure  2 for  y = 1.4.  Although  X(l)  = 4 
in  each  case,  the  functions  diverge  from  each  other  for  other  values  of 
m,  particularly  for  large  m,  and  the  detailed  shapes  of  the  functions  are 
rather  sensitive  to  y.  As  m tends  to  infinity,  the  three  functions  take 
on  limiting  values 


(4.20): 


A M = 2Y  ^ -1 

J y (Y  + l)  (y-D 


(4.28):  a2(")  = 3 TTf  (4.47) 

(4.34):  x3(»)  = 1 + f 

Clearly,  the  three  limiting  values  depend  strongly  on  y.  This  is  shown 
graphically  in  Figure  3.  Certainly,  the  behavior  of  s in  the  general  case 
in  (4.38)  must  also  depend  strongly  on  the  details  of  the  energy  function 
e(u,n) • 


1 / Mach 


Number  or  Amplitude 


Approximations 


The  preceding  discussion  has  been  based  on  the  asumption  that  one  of 
the  coupling  functions  is  known  exactly,  and,  if  that  is  the  case,  then  the 
development  of  the  equations  as  a hyperbolic  system  is  entirely  correct.  In 
an  actual  case,  the  coupling  functions  cannot  be  known  exactly  in  advance 
since  the  domain  of  dependence  for  the  full  equations  (3.3)  includes  a 
segment  of  the  shock  wave  itself.  This  occurs  because  the  backward  wave 
cone,  with  apex  on  the  rear  of  the  shock,  intersects  the  past  trajectory  of 
the  shock.  Thus,  the  hyperbolic  description  can  only  be  approximate  at 
best.  It  is  natural  to  ask,  therefore,  which  of  (4.20),  (4.28),  or  (4.34) 
will  provide  the  best  basis  for  an  approximate  treatment.  The  three  coup- 
ling terms  are  written  again  here  for  easy  comparison. 


(4.20): 

(4.28): 

(4.34): 


(4.48) 


In  each  case,  for  m = 1 , the  coupling  term  vanishes,  and  in  fact,  the 
governing  partial  differential  equations  all  become  the  same  ordinary 
differential  equation  along  a ray  because  the  curvature  terms  drop  out  as 
well.  As  the  shock  strength  increases  u decreases  from  1 in  the  range 
1 > F > U • > o,  a increases  from  o in  the  range  o < a < a <1,  and  l) 

increases  from  c in  the  range  c <.  U < 00 , as  may  be  seen  in  Figure  1 for 
the  special  case  of  a polytropic  gas  with  y = 1.4.  It  seems  clear  that  of 
the  three  coupling  terms  in  (4.48)  the  last  is  weakest  by  far.  For 

2 2 

example,  in  the  case  of  a polytropic  gas  the  term  (1-u  )/y  increases  mon- 

2 2 2 

otonically  to  (y + 1 ) / ( Y - 1 ) ( =6  for  y = 1.4);  the  term  m (l-a)(l-y  )/u 

2 

is  exactly  equal  to  m^-l  and  so  it  always  increases  towards  but  the 

term  (1-a) (l-y)/m  increases^ to  a weak  maximum  (about  0.11  at  m = 1.4  for 

y = 1.4),  then  decreases  to  zero.  A few  trial  calculations  indicate  that 

the  maximum  is  not  particularly  sensitive  to  and  never  exceeds  0.13  or 

so.  The  bracketed  terms  in  (4.48)  are  shown  in  Figure  4.  The  three  solid 

curves  correspond  to  the  case  y = 1.4.  Note  that  the  scale  for  C.  is 

expanded  to  forty  times  that  for  C.  and  C,,.  For  comparison  the  peak  value 

1 z 

of  C,  is  shown  on  the  same  scale  as  the  other  two,  and  appears  as  the 

5 2 

small  tick  mark  just  above  n = 0.6.  The  two  curves  with  dashed  segments 
represent  for  the  extreme  cases  y = 1 (lower  curve)  and  y = «>  (upper 
curve).  Only  the  end  points  and  the  maxima  were  actually  calculated  for 
the  latter  cases. 
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Figure  4:  Coupling  functions  vs.  p*'  for  y = 1.4.  Refer 

to  the  main  text  for  full  explanation. 


On  the  other  hand,  the  strength  of  the  coupling  term  depends  as  well 
2 - 1 

on  o . [ul,  Pc  [o  ],  or  (pc)  (p  , + pco  . u ,)  , and  an  estimate  of 

, II  |l  , t 

these  terms  should  be  made  in  any  application,  if  possible.  For  example, 
in  the  case  of  small  perturbation  of  a nearly  plane  shock,  the  third  com 
bination  of  terms  depends  only  on  overtaking  disturbances  (see  Reference  S 
p.  268  and  273).  Therefore,  a good  approximation  should  result  by  negleci 
ing  these  terms  in  cases  where  the  effects  of  inhomogeneities,  caustics, 
or  focusing  are  to  be  studied;  that  is,  in  cases  where  the  overtaking 
disturbance  is  absent  or  weak  compared  to  other  effects.  This  is  the 
approach  adopted  by  Whitham.  A good  approximation  should  also  result  n 
cases  where  the  overtaking  disturbance  can  be  estimated  by  linear  tin  "i . 
as  in  the  case  of  weak  shocks. 


2!) 


i 


Comparison  With  Whitham's  Shock  Dynamics 

In  a series  of  three  papers  (see  References  5,  6,  7,  and  8)  G.  B. 
Whitham  developed  a theory  of  shock  dynamics  that  was  intended  to  describe 
medium  to  strong  shocks  in  cases  dominated  by  geometry  or  short  time  local 
effects.  The  theory  was  based  on  two  purely  kinematic  relations,  and  an 
assumed  nonlinear  relation  between  Mach  number  and  the  area  of  a ray  tube. 
This  last  relation  was  derived  from  a one-dimensional  formulation  of  the 
perturbation  of  a shock  that  propagates  down  a channel  with  slowly  varying 
cross  section.  Equation  (4.36)  with  zero  right  hand  side  corresponds  to 
the  first  of  Whitham's  kinematic  relations.  The  second  may  be  found  by 
applying  Euler's  formula  for  the  derivative  of  a determinant  to  (2.7). 


Equation  (4.49)  is  equivalent  to  the  formula  (b./J)  . = o,  which,  if 

i , i 

specialized  to  the  case  of  normal  rays,  corresponds  to  the  form  originally 
used  by  Whitham.  In  any  case,  with  the  aid  of  (2.7),  (3.4),  and  (4.14), 
equation  (4.49)  may  be  reduced  to 


(A2) 


Un . 

i.i 


,3 


(6..  - n.n.) (u.  .) 
ij  i J ^ i,r 


(4.50) 


J.  K ) Z 

where  A = Jdet  ArAf  . (The  square  root  of  the  determinant  of  the  surface 

metric,  which  is  a measure  of  surface  area,  was  called  the  area  function 
by  Whitham  and  denoted  A.)  The  use  of  an  assumed  relation  between  the 

U 

area  function  and  the  shock  speed,  A = f(U),  in  equations  (4.36)  and 
(4.50),  both  with  zero  on  the  right  sides,  leads  to  a system  of  hyperbolic 
equations 

JL  I U I 

= o 


f-  u ,, 

rr^  U _ + U n. 

U ,3  1,1 


(4.51) 


n.  _ + A.ru  = o 
i,3  i ,r 


The  characteristic  speed  of  propagation  for  these  equations  is 

s _ 1 ( U f(U)  ) f m2-l  ) 

c+  V)-f,(u)| 


(4.52) 
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It  turns  out  that  X Cm)  in  (4.46)^  and  in  (4.52),  as  given  in  (8),  are 

identically  equal.  Thus,  the  full  equation  (4.34)  represents  a generaliza- 
tion of  Whitham's  shock  dynamics  to  fluids  that  are  nonconductors  of  heat 
and  that  have  generally  inhomogeneous  conditions  ahead  of  the  shock. 
Furthermore,  by  working  with  exact  equations,  it  has  been  possible  to 
retain  the  complete  coupling  term.  This  should  permit  error  estimates 
and  systematic  approximations  to  be  made. 
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APPENDIX 


The  inhomogeneous  terms  in  equations  (4.20),  (4.28),  and  (4.34)  are 
listed  here. 


(4.20): 
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